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Attempts to extend the capillary-wave theory of fluid interfacial fluctuations to microscopic wave¬ 
lengths, by introducing an effective wave-vector ( q ) dependent surface tension a e a(q), have encoun¬ 
tered difficulties. There is no consensus as to even the shape of a e g(q). By analysing a simple 
density functional model of the liquid-gas interface, we identify different schemes for separating mi¬ 
croscopic observables into background and interfacial contributions. In order for the backgrounds of 
the density-density correlation function and local structure factor to have a consistent and physically 
meaningful interpretation in terms of weighted bulk gas and liquid contributions, the background 
of the total structure factor must be characterised by a microscopic g-dependent length £(g) not 
identified previously. The necessity of including the q dependence of £(g) is illustrated explicitly 
in our model and has wider implications, i.e. in typical experimental and simulation studies, an 
indeterminacy in £( q ) will always be present, reminiscent of the cut-off used in capillary-wave theory. 
This leads inevitably to a large uncertainty in the q dependence of a e g(q). 

PACS numbers: 05.20.Jj, 68.03.Kn, 68.03.Cd 


Understanding the nature of the interface separat¬ 
ing coexisting fluid phases has provided fundamental 
insights into the properties of matter, including the 
necessity for attractive intermolecular forces. Of par¬ 
ticular importance is the role played by thermally 
induced interfacial wandering. Classical Capillary-Wave 
(CW) theory [1-7] provides a remarkably successful 
description of the fluctuation properties of fluid in¬ 
terfaces, leading directly to the concepts of interfacial 
roughness and entropic repulsion. These underpin 
our modern understanding of interfacial phase transi¬ 
tions. Recent extensions of CW theory that include 
a wavevector dependent tension attempt to describe 
interfacial fluctuations at microscopic scales. However, 
these have proved problematic [8-20]. For example, 
there is no consensus as to the sign of the rigidity 
coefficient characterising the proposed order q 2 correc¬ 
tion to the equilibrium surface tension. In the present 
paper, we provide a worked example, using a simple 
analytically solvable density functional theory (DFT), 
which illustrates that problems have arisen because 
previous analyses have failed to identify a q dependent 
microscopic lengtlrscale £(g) crucial for any treatment 
attempting to extend CW theory to shorter wavelengths. 


observables such as the parallel Fourier transform of the 
density-density correlation function G(z, z'\ q) and local 
structure factor S{z\q) = fdz' G(z, z'; q). Since, in the 
long-wavelength limit, fluctuations of f?(x) translate the 
density profile, CW theory implies that, in the interfacial 
region, 


G(z,z';q) 


P'(z)p'(z') 


crq z 


S(z; q) 


App'(z) 


crq z 


(i) 


for q <C A. Here, p(z) is the equilibrium density pro¬ 
file, A p = pi — p g is the difference in bulk coexist¬ 
ing densities, and we have set ksT = 1. Hence, for 
g <C A, we might anticipate that the total structure fac¬ 
tor S(q)= f^ l L dz S(z;q), behaves approximately as 

S(q) « LgSg + L[Si + (2) 


where L g and Li are the macroscopic sizes of the gas and 
liquid phases, which have (g-dependent) bulk structure 
factors S g and Si, respectively. In (2), the first two 
terms are the expected ’’background” arising from the 
bulk phases and the third is the ’’excess”, containing the 
Goldstone mode contribution. 


The central idea of CW theory is that the thermal 
excitations of long-wavelength undulations in the local 
height £(x) of a liquid-gas interface are resisted by the 
equilibrium surface tension a . Thus, in the absence of 
additional pinning effects, such as gravity, the thermal 
average of the fluctuations satisfies ( |i?(q )| 2 ) = ksT/aq 2 , 
where £(q) are the Fourier amplitudes of ^(x), with x = 
( x,y ). This result is valid only for q <C A, where A 
is a cut-off of order the inverse bulk correlation length. 
The small q singularity is also manifest in microscopic 


It is certainly natural to ask if the success of CW theory 
is limited only to long-wavelengths (gC A). A particular 
issue, hotly debated, is whether one can extend the the¬ 
ory to allow for a wavevector-dependent surface tension 
<r e ff(g) that replaces a in (1) and (2) and then meaning¬ 
fully apply this interfacial based description at the micro¬ 
scopic scale. For example, it is accepted that dispersion 
forces lead to a non-analytic term in the low q expansion 
that can be assimilated as <J e ff(q) = cr + Aq 2 hiq + ■ ■ ■ 
[12, 13, 20], similar to curvature corrections to the ten¬ 
sion of a spherical wall-fluid interface [21], Here, the 
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coefficient A > 0 is proportional to the coefficient of r~ 6 
in the interatomic pair potential, and is insensitive to 
the precise definition of the interface location. A more 
contentious question is whether or not one can identify 
a physically meaningful <r e s(q) over the whole range of 
wavevectors reaching to the inverse atomic scale. Two 
strategies are possible. The original suggestion was to 
identify £(x) from a given molecular configuration and 
then define a q-dependent tension from (|£(q)| 2 ) [12]. For 
this, one may now use sophisticated many-body defini¬ 
tions of the interface which go far beyond the notion of 
a local Gibbs dividing surface [16, 22-24], However, one 
cannot then infer directly the behaviour of the correla¬ 
tion function and structure factors because, away from 
the g<A limit, interfacial fluctuations no longer merely 
translate the profile [16, 20, 24], Moreover, one does not 
know what are the ’’bulk” contributions to measured ob¬ 
servables [16]. A different but related strategy, linked 
more directly to scattering experiments, is to start from 
the microscopic observables, such as the measured to¬ 
tal structure factor S(q), and identify an effective ten¬ 
sion from the excess contribution [17-19], i.e. one defines 
a eS (q) via: 


S ex (q) 


(A p ) 2 

CTeff(g)<7 2 ' 


(3) 


This approach has the advantage that a e g(q) does not 
depend on the definition of the interface position. How¬ 
ever, it still requires that we first write S(q) = S bg (q) + 
S ex (q ), and determine a suitable background S bg (q) in 
terms of S g and 5). In this paper, we use a simple density 
functional theory to illustrate that, if such a separation 
is imposed, one is required to write 


structure factors S g = 1 /6(K 2 +q 2 ) and Si = l/b(nf+q 2 ). 
Here, n g = \/ <\>"(pg)/b and m = \J4>"(pi)/b identify the 
bulk correlation lengths £ g = 1/k 9 and £; = 1 /m, which 
are implicitly temperature dependent. Minimization of 
H[p] leads to bp"(z) = <j>'(p), the solution of which deter¬ 
mines the free interfacial profile over the macroscopic in¬ 
terval [—L g , Li\. It is well known that the square-gradient 
theory does not incorporate the CW broadening of the 
density profile. However, the identification (3) only in¬ 
volves the difference in the bulk densities, and should be 
determined reliably even at mean-field level. The local 
structure factor then follows from the Ornstein-Zernike 
equation 

b +® 2 + l ft" s i z ;q) = 1 ( 6 ) 

We now consider a Double Parabola (DP) potential 
<t>{p) = 6k 2 (p - pg) 2 / 2 for p < po and <j>{p) = bnf(p - 
Pi) 2 / 2 for p > po, where p 0 = fp g + (1 - f)p t and 
/ = K g /(n g + Ki). The profile follows as 

. , ( Pg + (po - Pg)e~ K ^ zl if z < 0 

P ( z ) = < (7) 

{ Pi + ( Po - pi)e KlZ if Z>{) 

Here, we have located the origin such that p(0) = po, 
where the gradient p'{ 0) = A p/(£ g + £z) is largest. A 
straightforward calculation gives a = 6(Ap) 2 /2(£ g + £/). 
Henceforth, we set 6=1. The liquid/gas asymmetry 
means that the Gibbs dividing surface Zq , defined by 

Jdz ( p(z ) - pg - (pi - Pg)0(z - Z G )^j =0, (8) 


S b9 (q) = (L g + ((q)) S g + (L, - <(«)) S t (4) 

where £(g) is a microscopic q-dependent lengthscale 
which weights the bulk gas and liquid contributions. 
The properties of £(g) and the limitations placed on the 
accessibility of cr e g(q) from measured structure factors 
are described. We argue that, in general, an uncertainty 
in a e g(q) will always arise from that in £(g), implying 
that in experiments, simulations, and more realistic 
DFTs, it will be extremely difficult to identify a robust 
q-dependent tension from measurements of only p(z), 
and S(z-q) or S(q), beyond perhaps the leading-order 
terms of its small q expansion. 

We consider a square-gradient description of the inter¬ 
face based on the Grand Potential functional [4] 

H[p] = Jdr | ^ (Vp) 2 + <j>(p) | (5) 

The free-energy density <j>(p) is a double-well potential 
modeling the coexistence of gas and liquid phases be¬ 
low the critical temperature which have different bulk 


lies at Zq = £ g , away from the origin (see [25]). 9(z) 

is the Heaviside step function. The local structure factor 
follows from Eq. (6): 


S(z;q) 


S g + (5(0; q) - S g ) e - K s(«)l 2 l if 2 < 0 

St + iS^q)- Si)e~ K d<i)z if 2 > 0 


and attains its maximum value at the origin 


5(0; q) 


f _ Kg(q) Sg + Ki(q) Si _ 

q 2 (n g (q) + Kg)- 1 + (ni(q) + m)- 1 


( 10 ) 


where K g (q) 2 = k 2 + q 2 , and Ki(q) 2 = k 2 + q 2 . 
One may also determine G(z,z';q) [25]; for example, 
G(0,0;q) = S(0; q)/(K g (q)S g + Ki(q)Si). Note that, 
in contrast to the density profile, the decay of S(z; q) 
depends on q , meaning that the expression for S(z;q) is 
consistent with the CW expectation (1) only at small 
q , where K g (q) « K g and Ki(q) « kq Consequently, no 
single q-dependent surface tension cr e g(q) can replace a 
in all of the expressions (l)-(2). 
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We wish to repackage our results as G(z, z'\ q ) = 
G bg (z,z';q) + G ex {z,z'-,q) and S(z;q) = S bg {z-,q) + 
S ex (z-,q), where the separation is done consistently (for 
example, S bg (z-,q) = f dz'G b9 (z,z';q)), and is subject 
to a few necessary reasonable physical constraints [25]. 
Specifically, G bg (z,z';q) must be continuous, contain 
lengthscales only defined for the bulk fluids and decay 
to the appropriate liquid and gas correlation functions 
when 2 and z' are far from the interfacial region. These 
requirements ensure that S bg (z',q) is continuous which is 
also necessarry since, otherwise, S ex (z;q) would be dis¬ 
continuous, which would not match sensibly with (1). For 
example, it is not appropriate to suppose S bg {z\ q) jumps 
from S g to Si at z = Z G , say, since S ex (z;q ) would be 
inconsistent with CW theory even as q —> 0. This means 
that, within our model, the only allowed separation of the 
local structure factor, which is physically meaningful, is 
of the form 

f S„ + (S bg ( 0; q) - S g ) if 2 < 0 

S b9 {z;q) = \ 

{ Si + (S bg { 0; q)-Si) e~ Kl ^ z if z > 0 

(11) 

and S ex (z; q) = S ex (0; q) e~ K ^ z \ for z < 0, and sim¬ 
ilarly for z > 0 but with Ki(q). This way, the back¬ 
ground contribution varies continuously from S g to Si 
through the interface. Here, 5(0; g) = S bg (0;q) + 
S ex (0-,q) with consistency demanding that 5 &9 (0;g) = 
G bg (0, 0; q)(K g (q)~ 1 +/v;(g) _1 ) [25]. We note that, with¬ 
out loss of generality, one can always write 

G bg { 0, 0; q) = f G (q) G g + (1 - f G (q)) G l (12) 

where G g = 1/2 K g (q), Gi = 1/2 m(q) are the 2D Fourier 
transforms of the bulk gas and liquid correlation func¬ 
tions [20] and /g(<z) is the "fraction” of the bulk gas con¬ 
tributing to the background correlation function at the 
z = 0 plane. Similarly, we can always write 

S hs (0;g) = fs(q)S g + (l-fs{q))S, (13) 

where fs(q) is the ’’fraction” of the bulk gas contribut¬ 
ing to the background structure factor. The framework 
is completed by integration of S bg (z;q) over [—L g ,Li] 
determining that S bg (q) is of the form (4), where ([(g) 
is defined by analogy with the Gibbs dividing surface, 
Eq. (8): 

jdz( K S bg {z-q)-S g ~{S l -S g )6(z- C(g))) =0. (14) 

The properties of £(g) and a e g(q), obtained via (3), are 
linked as seen in the rigidity coefficient, defined from the 
low-g expansion <r e g(q) = a + K e g q 2 + ■ ■ ■ , which depends 
explicitly on £(0): 

K eff = -J- (^g ~ €g€i + - -^Z G C(O)^ (15) 


Note that for the case of Ising symmetry (£; = £ g ), the 
rigidity is uniquely determined as K e g = ba^/A [20] 
since £(g) = 0. However, with asymmetry, even the sign 
of the rigidity may be altered depending on £(0) and the 
asymmetry ratio K g /ni- 

We can now consider the merits and physical interpre¬ 
tation of different separation schemes by noting that fix¬ 
ing any one of fn(q), fs(q), or C(q) determines the other 
two. Consider setting f G (q) = 1/2, i.e. G bg (0,0-,q) = 
\{G g +Gi). This gives C(g) = |(K Z (g) _1 -K ff (g) _1 ). Alter¬ 
natively, setting f s (q) = 1/2, i.e. S bg (0;q) = ±(S g + Si), 
leads to ((q) = \{Ki{q)~ l — Kg(q)^ 1 ). These are not en¬ 
tirely implausible, but in both cases the respective ’’frac¬ 
tions” of bulk gas and liquid contributing to the back¬ 
ground correlation function and structure factor are com¬ 
pletely ad hoc. Thus, the choice f G (q) = 1/2 corresponds 
to / s = (3 K g (q) + Ki(q))/(A(K g {q) + K t {q))) which takes 
a markedly different value from 1/2 when K g K;. Sim¬ 
ilarly, fs(q) = 1/2 implies f G = m(q)/(n g (q) + m{q)), 
which is very different from 1/2 except for the case of 
Ising symmetry. The inconsistency can be avoided by 
requiring that fs(q) = /g(<?)) which determines fs(q) = 
K g{q)/( K g(q) + k i(q))- This is the only scheme for which 
the fractions of liquid and gas contributing to S bg (0;q) 
are identical to those for G b9 (0,0;<?) for all q, and leads 
to 

CO?) = KiO?)' 1 - Kg(q)- 1 . (16) 

Note that £(0) = Z G so that the weighting of S g and Si in 
S b9 (q) at q = 0 over the interval \—L gi L{\ is exactly the 
same as that of p g and pi in the total number of particles 
(per unit area) N = f dz p(z) = ( L g +Z G )p g +{Li-Z G )pi■ 
However, £(g) vanishes as q increases, which shifts the 
liquid-gas balance in (4) to z = 0, where the density 
gradient is maximum. The q dependence of the weight¬ 
ing lengthscale, which is also present for the other two 
schemes, has not been appreciated previously. Note that 
going beyond the present mean-field treatment by allow¬ 
ing for capillary-wave induced broadening of the inter¬ 
face has a minor effect: the result for Z G is unchanged, 
and the position of the maximum in the density gra¬ 
dient is only weakly shifted [25]. Within the present 
DP calculation, the result £(0) = Z G follows also when 
one notes that the condition fs(q) = fc(q) implies that 
S bg (z;q) is continuous and differentiable at the origin. 
Thus, S bg (z; 0) has exactly the same shape as the den¬ 
sity profile. This also means that all other separation 
schemes lead to a S b9 (z\q) which has a kink at the ori¬ 
gin. Results for different choices of £(q) are shown in 
Fig. 1 of [25]. The effective tension follows as 

Veff(q) = 2 K g (q) 2 Ki(q) 2 / 1 + 1 \ 

<7 K g Ki(Kg(q) + Ki(q)) \K g (q) + Kg Ki(q')+KiJ' 

(17) 
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FIG. 1. Variation of the effective surface tension a e g(q) for 
increasing liquid/gas asymmetry obtained for different choices 
of the microscopic background weighting lengthscale £( q ) = 
e(Ki(g)“ 1 —K g (g) _1 ), and for £( q ) = Zg■ In the case of Ising 
symmetry (top left) all results are identical. According to 
simple DFT estimates [26], K g = 2 m corresponds to T/T c « 
0.75, and K g = Am to T/T c « 0.6. 


yielding a rigidity K eS = (3£ 2 - £ g £ t + 3£f) cr/4, which 
remains positive for all values of k 9 /k ;. 

Next, consider schemes which instead set £(g) to be a 
(^-independent) constant. Consider, for example, setting 
£(g) = Zq. While this initially appears desirable, it re¬ 
quires that S b9 (0;q ) is lower than S g (even negative) for 
sufficiently large q (see [25]). This arises as the reasonable 
requirement S g < S b9 (0;q) < Si leads naturally to the 
bounds -K g (g) _1 < ((g) < Kj(g) -1 , implying that ((q) 
must vanish as q increases. A final choice £(g) =0 corre¬ 
sponds to the plane where S(z;q) (and p'(z)) is largest. 
In this case, there is no problem with the behaviour of 
S b9 (0;q), for which f s (q) = Ki(q)/(K g {q) + K t (q)). 
However, this choice is equivalent to setting 
fo(q) = (2 Ki{q) - Kg{q))/(K g (q) + m{q)) in (12) which is 
negative for even moderate liquid/gas asymmetry. Thus, 
when K g ki, the fraction of the bulk gas contributing 
to the background correlation function is /g(<?) ~ — 1 
while the fraction for the liquid is 1 — /c(g) ~ 2. In 
Fig. 1, we show plots of a e g(q) obtained for different 
separation schemes and for increasing asymmetries 
Kg/m. For the case of pure Ising symmetry (k 9 = «/), 
all are equivalent and ([(g) = Zq = 0. When K g > ki , 
the difference with (17) is always largest for the choice 
£( q ) = Zq , which underestimates significantly cr e g(q) at 
large q, even though it identifies the rigidity correctly. 
The opposite is true for the other separation schemes 
which work reasonably well at large q but, alas, have 
different rigidities, consistent with Eq. (15). 


In summary, explicit DFT results for the simple DP 
model, tell us that while in principle there is freedom to 
choose any separation, S(q) = S b9 (q) + S ex (q), in all bar 
one case there was no physical interpretation of what the 
’’background” means in terms of weighted bulk contri¬ 
butions. In the ’’physical” case, Eq. (16), the weighting 
lengthscale £(g) is g-dependent. Let us consider the wider 
implication. Suppose that for another model, simulation 
or experimental study we have two different separation 
schemes with different £( q ) and a e g(q). Since the total 
structure factor S(q) is the same in both descriptions, 
the difference in the inverse g-dependent surface tensions 
follows from (3) and (4): 

A (^faj) = ^ (Pi- PgY AC(9) (18) 

where A£(g) is the difference in weighting lengthscales. 
The two schemes agree the value of a c g(q) in the limit 
q —> 0 must be a, but disagree at larger q, though the 
difference A£)(g) remains microscopic. The expression 
(18) can therefore be viewed as characterising a funda¬ 
mental uncertainty in the g-dependent surface tension 
arising from the indeterminacy of C,(q). This would not 
be a problem if the weighting lengthscale ((q) was equal 
to the position of the Gibbs dividing surface for all q. 
However, the fact that £(g) is g-dependent, even in this 
very simple DFT, means that this must be the case 
more generally. Therefore, unless ((q) can be measured 
or determined independently , there will always be an 
uncertainty in the excess contribution to the structure 
factor implying that the behaviour of cr e ff(g) away from 
g = 0 is essentially unknowable. In this regard, the 
limits of certainty on the form of cr e g(q) arising from 
an unknown microscopic lengthscale are reminiscent 
of the high wavevector cutoff A used in classical CW 
theory. One arrives at very similar conclusions if one 
defines a g-dependent tension via the local structure 
factor, instead of (3). For example, one could measure 
where S(z;q) is largest (in our case z = 0), and define, 
S es (0; g) = App l (0)/q 2 a e g(q) which generalises (1). 
In this case, we are still left with an uncertainty in 
creff(g), similar to (18) except that A£(g) is replaced by 
(£ g + £;)A/s(g), where A/s(g) is the uncertainty in the 
g dependence of the weighting fraction /s(g). 

In our analysis, we found that £(0) = Zq-, this result 
was not imposed. Rather, it emerged from trying to 
identify a consistent choice for the fractions /s(g) and 
/ G (g) in the background S bg (0;q) and G b9 (0,0;g). It 
would be extraordinary if, beyond the present DP 
model, a suitable separation of the local structure factor 
always results in a weighting lengthscale satisfying 
£(0) = Zq. One way of imposing this would be to 
set S b9 (z; 0) - S g oc p{z) — p g . Although this is valid 
for the separation leading to (16) in the present DP 
model, it cannot be generally valid. Recall that, beyond 
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mean-field, the profile is strongly affected by thermal 
wandering which cannot, by definition, be in the back¬ 
ground contribution. This leaves us with two scenarios 
which determine the robustness of the q expansion of 
a e ff(q). First, if as found here, sensible separations 
always find £(0) = Zq, then the uncertainty A ((q) van¬ 
ishes as q —> 0. From (18), this means that one may then 
write cr e ff(g) = a + A q 1 2 3 4 5 In q + K e s q 2 + • • • and identify 
a meaningful rigidity K e g. Nothing else, however, can 
be said unless the q dependence of £( q ) is determined. 
Second, if one finds instead that there are, in general, 
different equally acceptable ways of separating S(z;q) 
into background and excess terms which sometimes 
result in £(0) Zq, then the uncertainty A<((0) yL 0 
implies that only the term 0(q 2 lng) induced by the dis¬ 
persion forces remains well characterised. In this second 
scenario, which appears more likely to us, it is not just 
that AT e ff is non-unique but that the separation of S(q) 
into background and excess terms may well be ill-defined. 
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Here we provide further details of three calculations: 
The determination of G(z, z '; q), its separation into back¬ 
ground and excess terms, and the influence of capillary- 
waves fluctuations on lengthscales characterising the in¬ 
terfacial profile. A plot showing the background local 
structure factor S bg (z; q) for different choices of ((q), and 
the location of the Gibbs dividing surface Zq is also in¬ 
cluded. 


I. THE EXPRESSION FOR G(z,z'-,q) 

G(z , z'; q) is obtained from solving the Ornstein- 
Zernike equation which, for our local square-gradient 
functional, reduces to [1] 

(~ +9 2 +^"(p(^))) G{z,z'-,q) = S(z - z') (1) 

where, as earlier, we have set ksT = 1 and 6=1. Far 
from the interfacial region, the correlation function tends 
to the appropriate bulk expression: 

-K a (g) \z-z '| 

< 2 > 

on the gas (a = g) and liquid sides (a = l ) of the interface, 
respectively. Within the DP model, the second derivative 
of the potential can be written <t>"{p(z )) = k 2 + (/t 2 — 
Ac 2 ) Q(z) — (n g + Ki) S(z — z'), and the Ornstein-Zernike 
equation is readily solved using elementary methods. The 
correlation function takes its maximum value when both 
particles are at the origin (which have located at the point 
where the derivative p'(z) is maximal), and is given by 

G(0,0;5) = ^—:;-—r-r~r (3) 

q 2 {n g (q) + K g ) 1 + (m{q) + Hi) 1 

Thus, the correlation function exhibits the anticipated 
Goldstone mode singularity as q —> 0. Then, if the par¬ 
ticles are on the same side of the interface ( zz 1 > 0), 
G{z, z'\ q) is given by 

G = G a (\z - z'\\ q) + (G(0,0; q) - G a ( 0; g)) e -".(*)(l*+* , l) 

(4) 

with a = g or a = l, as appropriate. However, if these lie 
on either side (say, z < 0 < z'), 

G = G(0,0;q)e- K » (9)| d e -«d9)P (5) 

The correlation function is a continuous function of z and 



FIG. 1. Background local structure factor for three differ¬ 
ent choices of the microscopic weighting lengthscale (~(q ) = 
e(«((g)“ 1 — K g (q')~ 1 ), and for ((q) = Zg (red line). Here, 
Kg/m = 2, and q = 2 k;. For e ^ 1, corresponding to schemes 
where fs fa, the S bs (q; z) has a kink at the origin. For the 
scheme with ((q) = Zg, S b9 (z;q) is lower than S g for z < 0. 
The inset shows the density profile and the location of the 
Gibbs dividing surface, Zg- 


II. THE SEPARATION OF G(z,z'-,q) 
A. Ising symmetry: m = K g 


The expression for G(z, z'q) simplifies considerably 
when ki = K g = k (i.e. Ising-like symmetry). The bulk 
liquid and gas phases have identical correlation functions 
and structure factors which we write as 


G b {\z-z'[,q) 


e -K(q)\z-z'\ 

2s{q) 


S\q) 



( 6 ) 


where n(q) = \Jk 2 + q 2 . In this case, the position de¬ 
pendence of the correlation function throughout the in¬ 
homogeneous region can be expressed as [2] 


G(z, z'\ q) = G b (\z - z'\;q ) 


G b (\z\-g)G b (\z'\;q) 
G b (0; 0) — G b (0; q) 1 ’ 


This result splits unambigously into a background term, 
equal to the bulk correlation function, and an excess part 


_ G b {\z\-q)G b {\z'\-q) 
1 ’ ,q) G b {0:0)-G b {0;q) 


(8) 
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which decays exponentially, as each particle moves away 
from the interface, and is controlled by the inverse length- 
scale n(q). For q <C k, and distances close to the interface 
(8) is entirely in keeping with the prediction of capillary- 
wave theory (Eq. (1) in the main article). It follows that, 
for this Ising symmetric case, the local structure factor 
also separates [2] 


, . S l (9)C l (| 2 |;9) 

S( "' 5) S (?) + G b (0;0) — G' b (0;</) (9) 


so that the background contribution is simply equal to 
the bulk Sb(q). Similar to the correlation function, the 
excess contribution decays exponentially as the particle 
position 2 moves away from the interface. Integration of 
S(z;q) over the macrosopic interval [— L g ,L{\ gives 


Stot(q) 


( L g + Li) S b (q) + 


S b (q) 2 

G b (0; 0) — G h (0; q) 


( 10 ) 


if zz' > 0 and 

Qhg = G b 9 ( 0i 0 . e ~Kg(.q)\z\ e -Ki(q)z' ( 14 ) 

if 3 < 0 < z'. As for the symmetric Ising case (see (8)), 
the excess term is a simple product 

G ex (z, z'; q) = G ex { 0,0; q) e ~ K ^ z '^ (15) 

This decomposition reproduces identically the equi¬ 
librium function G(z,z;q ), provided that we impose 
G(0,0; q) = G b9 (0,0;g) + G ex { 0,0; g). This leaves us 
with only a single function characterising the whole 
separation: the value of G bs (0,0;g). Integration leads 
directly to the separation of the local structure factor- 
described in the main article and to the consistency 
requirement S bg (0;q) = G bg (0, 0; q)(K g (q)~ 1 + Ki(q)^ 1 ). 
Thus, one may equally consider S b9 (0;q) as the function 
to characterise the separation; see Eqs. (12) and (13) of 
the main article. 


where the first term is the total background contribution. 
From Eq. (4) of the main article, it follows that ((q) = 
0, and one may then identify a wave-vector dependent 
tension from the excess contribution as 


(Jeff (q) 


2(1 + g 2 4 2 ) 2 

1 + + V 71 + 


( 11 ) 


where £& = 1/k is the bulk correlation length. This is the 
result given in Eq. (28) of [2]. 


We emphasize that any separation/decomposition 
of the local structure factor S(z; q) into background 
and excess, which dictates the separation of the total 
structure factor S(q) as defined in the main article, 
must be consistent with the corresponding separation of 
G(z, z’\q). It is the latter, the density-density correlation 
function (clearly, the fundamental statistical mechanical 
quantity) which describes the structure of the interface 
at the two-particle level. 


B. Liquid-gas asymmetry: ki ^ K g 

We wish to see if such a simple separation of G(z, z; q) 
occurs also when there is a more general asymmetry be¬ 
tween the bulk phases. That is, we wish to write 

G(z, z 1 ; q) = G b9 (z, z;' q) + G ex (z , zf q) (12) 

where now the background ( bg ) contribution necessarily 
involves a mix of bulk liquid and gas correlation func¬ 
tions. It is only permissible to construct the background 
and excess contributions using the same elementary expo¬ 
nential functions e ±K “G) z and appearing in the 

equilibrium solution (4) and (5). Otherwise, one is intro¬ 
ducing entirely arbitrary and non-physical lengthscales 
which are neither incorporated in the bulk nor in the in¬ 
homogeneous system. It is also necessary that the back¬ 
ground and excess contributions are continuous functions 
of z and z'\ otherwise, there is no connection with the 
standard capillary-wave expression (Eq. (1) in the main 
article) as q —> 0 nor with effective Hamiltonian theory. 
Finally, the background must approach the appropriate 
bulk function far from the interface. This leaves us only 
with a background contribution G b9 (z,z';q) of the form 

G b9 = G a (\z-z'\-,q)+(G b9 (0, 0; g)-G o (0; q)) e-^XM+l-'D 

(13) 


III. THE INFLUENCE OF CAPILLARY-WAVE 
FLUCTUATIONS ON THE DENSITY PROFILE 

Beyond mean-field approximation, the equilibrium 
profile Peq(z) is broadened by interfacial fluctuations. 
This is well-described by the approximate capillary-wave 
formula [3] 

Peq(z) « jdl p 0 (z - i) P{tj (16) 

where here po(z) may be taken to be the underlying 
mean-held profile (Eq. (7) in the main article) and P(£) 
is the probability distribution for the interfacial position 
i. If the interface has finite transverse area Ljj, and there 
are no other pinning fields, then P{€) has a simple Gaus¬ 
sian form P{£) = e~ t ‘ where the interfacial 

width £j_ satisfies the celebrated capillary-wave formula 
£)[ = (27ro-) -1 ln(L||A) (recall that cr is the surface ten¬ 
sion, and A the cut-off). In the thermodynamic limit, 
L || —► oo, the equilibrium profile is very different from 
the underlying mean-held prohle, and is well-described 
an error function of width £j_ [3]. However, the posi¬ 
tion of the Gibbs dividing surface Zq and the location 
of the maximum in the density gradient Zq are hardly 
affected. In fact, the normalization condition on P(£) 



means that the location of Zq is completely unaltered and, therefore, remains close to the origin as Ly —> oo 
from its mean-field prediction, Zq = £,1 — £ g . Similarly, similar to the mean-field prediction, 
asymptotic analysis of the integral shows that the max¬ 
imum in the density profile (which occurs at Zq = 0 in 
mean-field approximation) now occurs at 




Zcitt + Zf) 

a 


(17) 
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